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(57) Abstract 



In the method for constructing a fault dictio- 
nary, a description of the circuit is first made. A list 
of faults is extracted from this description of the cir- 
cuit, and a fault-free circuit distribution of an output 
parameter of the circuit is calculated in response to 
the circuit description. A faulty circuit distribution of 
the output parameter is calculated in response to the 
faults of the list, and a fault value is calculated from 
the fault-free and faulty circuit distributions. The cal- 
culated fault value is stored in the fault dictionary in 
view of subsequently specifying a test vector for ap- 
plication to the circuit in view of testing this circuit. 
To specify test vectors in view of testing a given cir- 
cuit, (a) a set of stimuli is first selected and then (b) 
a stimulus of this set is selected, (c) a fault from a 
list of faults of the circuit is selected, (d) a fault value 
related to the selected fault and stimulus is found in a 
fault dictionary, (e) the fault value found in step (d) is 
compared to a typical fault value, (f) whether the se- 
lected fault is detected or undetected by applying the 
selected stimulus to the circuit including the selected 
fault is determined in accordance with the result of 
the comparison of step (e), (g) steps (b, c, d, e and f) 
are repeated for each pair of stimulus and fault; and 
(h) test vectors are built from the stimuli of the set 
most susceptible to detect faults in the circuit. 
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METHOD FOR PARALLEL ANALOG AND DIGITAL CIRCUIT 
FAULT SIMULATION AND TEST SET SPECIFICATION 

5 

BACKGROUND OF THE INVENTION 
1. Field of the invention: 

10 

The present invention relates to a highly effective 
method for conducting parallel fault simulation and test set specification 
of digital and analog circuits. The method according to the invention 
formulates the fault simulation problem as a problem of estimating the 
15 fault value based on the distance between fault-free and faulty circuit 
distributions of the butput parameter. 



2. Brief description of the prior art: 

20 

Previous work in fault simulation and test generation 
focuses on digital circuits using the classical stuck-at fault model. Both 
serial and parallel fault simulation techniques have been developed. 
Algorithms for serial fault simulation constitute the simplest method of 
25 simulating faults. It consists of transforming the model of a fault-free 
circuit N to model a circuit Nf created by the fault f. Then Nf is simulated. 
The entire process is repeated for each fault of interest. Thus faults are 
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simulated one at a time [M. Abramovici, M.A. Breuer, A.D. Friedman, 
"Digital Systems Testing and Testable Design", IEEE press, 1990]. Other 
fault simulation techniques, for example parallel, deductive, and 
concurrent fault simulation techniques have been developed and differ 
from the serial method in two fundamental aspects: a) they determine the 
5 behaviour of the circuit N in the presence of a fault without explicitly 
changing the model of the fault-free circuit N and, b) they are capable of 
simultaneously simulating a set of faults. 

Research in the area of analog circuit fault simulation 
10 and test vector generation has not been as successful as its digital 
counterpart owing to the difficulty in modelling analog behaviour, the 
continuous nature of the analog input and output signals, the non-linearity 
of the circuit elements, and the complicated relation between the input 
and output signals referred to as "transfer function". Thus, because of the 
15 complex electrical nature of analog circuits, a direct application of the 
digital models has proved to be inadequate for capturing the faulty 
behaviour. 

In the article of L. Milor and V. Visvanathan entitled 
20 "Detection of Catastrophic Faults in Analog Integrated Circuits", IEEE 
Trans. Computer-Aided Design, Vol. CAD-8, no.2, pp. 114-130, Feb. 
1989, it has been suggested to model the faulty analog behaviour as a 
modification to the nominal macromodel. For instance, the fault model for 
a transistor has been implemented by replacing each transistor by a 
25 transistor surrounded by switches as shown in Figure 1. A faulty circuit 
can be obtained from the fault-free circuit by opening or closing the 
appropriate switch(es) 10,11,12 and/or 13. 
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In the article of R.J.A. Harvey, A.M.D. Richardson, 
E.M.J.G. Bruls, K. Baker, entitled "Analog Fault Simulation Based on 
Layout Dependent Fault Models", ITC 95, pp-641-649, to enable the 
circuit fault-effects to be simulated with a reasonable simulation time, 
behavioural models for each circuit block were developed. Hybrid fault 
5 simulations were performed by replacing each circuit block with its 
equivalent behavioural model to insert a target fault. Each fault is 
manually inserted in a netlist for simulation. The behavioural models 
reduced the simulation time by a factor of 10 to 36. 

10 Hard fault modelling and simulation has been the subject 

of many publications: 

- R.J.A. Harvey, A.M.D. Richardson, E.M.J.G. Bruls, K. Baker, "Analog 
Fault Simulation Based on Layout Dependent Fault Models", ITC 95, pp- 

15 641-649; 

- L. Milor and V. Visvanathan, "Detection of Catastrophic Faults in Analog 
Integrated Circuits", IEEE Trans. Computer-Aided Design, Vol. CAD-8, 
no.2, pp. 114-130, Feb. 1989; 

20 

- Naveena Nagi, A. Chatterjee, Jacob A. Abraham, "Fault Simulation Of 
Linear Analog Circuits", Analog Integrated Circuits and Signal Processing 
1993; 

25 - Majoj Sachdev, "A Realistic Defect Oriented Testability Methodology for 
Analog Circuits", JETTA 1993; and 
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- Naveena Nagi and Jacob A. Abraham, "Hierarchical Fault Modeling For 
Analog and Mixed-Signal Circuits", IEEE VLSI Test Symposium 1992, pp 
92-101. 

The presented approaches are all based on cause-effect analysis and do 
5 not enable parallel fault simulation. Indeed, cause-effect analysis 
enumerates all the possible faults (causes) existing in a fault model and 
determines all their corresponding responses (effects) to a given applied 
test in a serial manner. The required simulation time can become 
impractically long, especially for large analog designs. 

10 

OBJECTS OF THE INVENTION 

15 An object of the present invention is to provide a method 

for performing parallel fault simulation and test vector specification based 
on effect-cause analysis. From the distance between a fault-free circuit 
distribution and a faulty circuit distribution (effect), the method according 
to the present invention approximate a fault (cause) value for all modelled 

20 faults simultaneously by linear estimation. 

SUMMARY OF THE INVENTION 

25 

More specifically, in accordance with the present 
invention, there is provided a method for constructing a fault dictionary for 
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a circuit in view of subsequently testing the circuit, comprising the steps 
of: 

( a) describing said circuit; 

( b) extracting a list of faults from said description of the 

circuit; 

5 ( c) in response to the description of said circuit, 

calculating a fault-free circuit distribution of an output parameter of the 
circuit; 

( d) calculating a faulty circuit distribution of the output 
parameter in response to the faults of said list; 

10 ( ©) calculating a fault value from the fault-free and faulty 

circuit distributions; 

( f) storing the calculated fault value in the fault 
dictionary in view of subsequently specifying at least one test vector for 
application to the circuit in view of testing said circuit. 

15 

The present invention also relates to a method for 
building test vectors in view of testing a given circuit, comprising the steps 
of: 

( a) selecting a set of stimuli; 
20 ( b) selecting a stimulus of said set; 

( c) selecting a fault from a list of faults of said circuit; 

( d) from a fault dictionary, finding a fault value related 
to the selected fault and stimulus; 

( e) comparing the faut value found in step ( d) to a 
25 typical fault value; 
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(f) determining whether the selected fault is detected or 
undetected by applying said stimulus to the circuit including the selected 
fault in accordance with the result of the comparison of step ( e); 

( g) repeating steps ( b), ( c), ( d), ( e) and ( f) for each 
pair of stimulus and fault; and 

5 ( 9) building test vectors from the stimulus of said set 

most susceptible to detect faults in said circuit. 

The objects, advantages and other features of the 
present invention will become more apparent upon reading of the 
10 following non restrictive description of a preferred embodiment thereof, 
given by way of example only with reference to the accompanying 
drawings. 

15 BRIEF DESCRIPTION OF THE DRAWINGS 

In the appended drawings: 

20 Figure 1 , which is labelled as "prior art" is a fault model 

for a transistor, consisting of a transistor surrounded by switches; 

Figure 2 is a schematic view of an inverter layout, 
including an intermediate list of faults of this inverter layout, and a final list 
25 of faults obtained from the intermediate list after redundancies have been 
eliminated; 
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Figure 3 illustrates the circuit of a simple voltage divider; 

Figure 4 is a complete list of 9 faults that can be found 
in the voltage divider of Figure 3; 

5 Figure s is a reduced list of 6 faults built from the 

complete list of Figure 4 after removal of redundant and undetected 
faults; 

Figure 6 is a graph showing examples of fault-free and 
10 faulty circuit distributions; 

Figure 7 is a fault simulation flow chart as implemented 
by the method according to the invention for constructing a fault 
dictionary; 

15 

Figure 8 is a hard fault test vector specification 

algorithm; 

Figure 9 is a schematic diagram of a second order band- 

20 pass filter; 

Figure 10 is a list of 5 faults connected to the second 
order band-pass filter of Figure 9, indicating that 4 out of the 5 faults have 
been detected by means of a single test vector; 

25 

Figure 11 is a graph of the output parameter of the 
second order band-pass filter of Figure 9 versus the frequency of the test 
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vector, showing that for a test vector of 10 kHz the 4 detected faults of 
Figure 10 modified the output voltage (output parameter) by more than 
5% with respect to the fault-free output voltage; 

Figure 12 is a schematic diagram of the circuit of a fifth 
5 order chebychev filter; and 

Figure 13 is a graph showing the fault coverage as a 
function of the number of test vectors for an inverter, a low pass filter, a 
state variable filter, a chebychev filter, a 4-bit analog-to-digital converter, 
10 an 8-bit analog-to-digital converter flash, and an 8-bit current digital-to- 
analog converter. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

15 

Fault simulation is used to construct a fault dictionary. 

Conceptually, a fault dictionary stores the signatures of faults in response 

to a specific stimulus T. Prior art fault simulation techniques are fault 
20 oriented and consists of: 

( a) transforming a fault-free circuit N to a faulty circuit Nf 

created by a fault /to simulate Nf (the fault required by 
the fault-free circuit N to create the faulty circuit Nf could 
be obtained from statistical analysis of circuit defects 
25 (Table 1)); 

( b) simulating the circuit in the presence of the fault f, and 
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( c) storing the signature (value of the output parameter) of 

each fault in the fault dictionary for use in test vector 
specification. 

Constructing a fault dictionary through computation of every possible fault 
5 before testing is impractical. 



TABLE 1 Upper and lower resistances used for hard fault modelling 



10 





Lower 


Upper 




Resistance 


Resistance 


Defect type 


(ohms) 


(ohms) 


Added metal 1 


0.2 


1000 


Added metal 2 


0.2 


1000 


Via short 


5 


5 


Junction leakage 


100 


10 000 


Poly-metal 1 short 


0.2 


1000 


Poly-metal 2 short 


0.2 


1000 


Poly-poly short 


20 


1000 


Open 


1 Meg 


oo 



20 

Figure 7 is a fault simulation flow chart, as implemented 
in the method according to the invention. 

Step 1 9: A layout and/or schematic description of the circuit to be 

25 tested is first made. 
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Step 20: From the description of step 1 9, a layout-based fault list 

and/or a schematic-based fault list is/are extracted to 
form a fault list 21. 

Step 22: The sensitivities of the output parameter to variations of 

5 circuit components due to manufacturing process 

variation is computed from the circuit description of step 
1 9. The sensitivity is defined as the ratio of the relative 
deviation of the output parameter to the relative 
deviation of the circuit component. 

10 

Step 23: In step 23, the fault-free circuit distribution (see Figure 

6) of the output parameter is computed in response to 
the circuit description (step 19) and the computed 
sensitivity (step 22). 

15 

Step 24: A stimulus T is selected. 

Steps 25-26: The gradients of the output parameter with respect to all 
faults (step 25) are computed for the selected stimulus 
20 T (step 24). The output of step 26 is representative of 

the faulty circuit distribution (see Figure 6) of the output 
parameter in response to stimulus T. 

Step 27: The fault value R fault is calculated from the gradients 

25 computed in step 26 and the fault-free circuit 

distribution computed in step 23. 
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Step 28: calculated in step 27 is stored in the fault dictionary 

to form this fault dictionary. 

Step 29: If the response to additional stimuli is desired, the 

algorithm returns to step 24. 

5 

Step 30: If no additional stimulus response is desired (step 29) 

the algorithm terminates. 

0 Several steps of the algorithm of Figure 7 will be more 

clearly described hereinafter. 

In accordance with the present invention, the fault 
dictionary is not constructed by storing the output signature of the fault f 
(effect), but by computing and storing the value R fautt of the fault (cause) 
by parallel fault simulation. For instance, R^, indicates the value of the 
resistance that, if added to the circuit, will drive the output parameter 
under test out of a predetermined tolerance range. The output parameter 
is defined as a measured performance; the output parameter may be a 
voltage amplitude, a current amplitude or any other circuit response or 
specification. The fault list is first extracted and the fault value R fault 
associated to the selected stimulus is then simultaneously calculated for 
all faults of the list from the fault-free and faulty circuit distributions of the 
output parameter, and stored in the fault dictionary. As will be explained 
in the following description, fault coverage and test vector specification 
can also be performed. 
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Steps 20 and 21: 

The method according to the invention use a layout- 
based fault list and/or a schematic-based fault list. The layout-based fault 
list is used at macro level while the schematic-based fault list is used at 
5 a higheMeveli Le. at the level of interconnections between modules of the 
circuit. The schematic-based fault list is also used at early design stage 
where macro layouts are not available. 

The method used to generate the layout-based fault list 
1 0 consists of moving a defect over the entire area of the layout. At each 
position taken by the defect, the polygons touch and their net numbers 
are obtained to deduce the fault produced by the defect [H. Wlaker and 
S.W. Director, "VLASIC: A Catastrophic Fault Yield Simulator for 
Integrated Circuits", IEEE Transactions on Computer Aided Design of 
15 Integrated Circuits and Systems", Vol. CAD5(4), pp 541-556, Oct. 1 986]. 
When the defect has been moved over the entire layout, the layout-based 
fault list is complete. Each fault is described by its type , the coordinates 
of the defect responsible for the fault, and a probability related to the size 
of the responsible defect. 

20 

Example: A fault list was obtained for the inverter layout 14 of 

Figure 2. A total of 240 defects were inserted in the 
layout 14, which resulted in 26 faults. An intermediate 
faut list 15 gives the type of the fault and the 
25 coordinates (position) of the defect causing the fault. 

Eliminating the redundancy (step 17) led to the final 
fault list 16 of Figure 2, which contains three different 



WO 98/55880 




PCT/CA98/00538 



13 



faults. In the final fault list 16, each fault is weighted by 
the number of occurrences. 

« 

Extraction of a schematic-based fault list is illustrated in 
Figures 3-5. The list shown in Figure 4 is not exhaustive since the only 
5 shorts-listed are between the same element nodes. The method 
according to the invention also enables generation of all combinations of 
two, three or more shorted nodes. Using this feature may lead to a large 
list of faults containing a significant amount of nonrealistic faults. The 
example treated in Figures 3-5 does not include these combinations. 

In Figures 3-5, a simple voltage divider 18 (Figure 3) is 
taken as an example. The complete list of Figure 4 
contains 9 faults. After having withdrawn the 
undetected and redundant faults from the complete list 
of Figure 4, the reduced list of Figure 5 containing 6 
faults is obtained. 

As an additional example, for a MOSFET transistor, the 
fault list is constructed of three shorts and two opens: short gate-source, 
short gate-drain, short source-drain, open drain, and open source. 

Fault modelling and fault dictionary generation: 

As indicated in the foregoing description, the present 
invention proposes to construct a fault dictionary by storing the fault value 
RfcuH (cause) instead of storing the output parameter (effect). The effect 
is a constant value which represents the detectability threshold. The 



15 



Example: 
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cause is the minimal fault value that, if added to the circuit, will cause the 
output parameter to go out of a predetermined tolerance range and make 
the fault detectable. One way is to define the detectability threshold as 
( i) a constant absolute value deviation of the output parameter from its 
nominal value, for example D thKshM = 10 mv, ( ii) a constant percentage 
5 of the fault-free value of the output parameter, for example D tmshold = 5%, 
or (iii) a constant factor of the fault-free circuit distribution of the output 
parameter and a constant factor of the faulty circuit distribution of the 
same output parameter, for example 3o ff and 3o>, respectively. 

The first two methods ( i) and ( ii) to define the 
detectability threshold are straightforward. 

Steps 22 and 23: 

For the third method ( iii) (constant factors of the fault- 
free circuit distribution), a piecewise linear estimation is used to compute 
the fault-free circuit distribution of the output parameter due to process 
variations using Equation (1) [Nairn B-Hamida, Khaled Saab, David 
Marche and Bozena Kaminska, "FaultMaxx: A Perturbation Based Fault 
Modeling and Simulation for Mixed-Signal Circuits", IEEE Asien Test 
Conference, October 1997; and Nairn B. Hamida, Khaled Saab, David 
Marche, B. Kaminska and Gay Quesnel, "LIMSoft Automated Tool for 
Design and Test Integration of Analog Circuits", International Test 
Conference, 1996]. From Equation(1), the mean and standard deviation 
of the fault-free circuit output parameter can be obtained by means of 
Equations (2) and (3), respectively. 
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out = out, * £ S * Ax, (1) 

M 



Mo* = Pout. + E S x ° u ' MaX( (2) 

/=1 



^-EKWx.-EEsC's-o (3) 

/=1 /*yM/=l ' ' 



where - out 0 is the nominal value of the output parameter; 

- out is the estimated value of the output parameter due 
5 to variation of the components; 

-Ax, is the variation of the circuit component x, due to 
process variation; 

- S xl out is the sensitivity of out to x,; 

- o,,, is the standard deviation of the component x,; and 
10 - o^- are covariance terms. 

Steps 24-26: 

Now, if the faulty circuit distribution of the output 
15 parameter and its standard deviation due to added resistance are 
computed, equations (1), (2) and (3) respectively become 
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— IW---M^ t s:; m AX( + - 6,;«" (5) 



2 



where R /; is the newly added component due to short or open and G ff f* 
is the gradient of out with respect to the fault /, , is the standard 
deviation of the fault value (resistance), and o^,, is the covariance term 
5 between the newly added component and the components in the original 
fault-free circuit. 

Since o Rfi is always equal to zero, and if we consider the 
variable as independent, the covariance terms o^ are equal to zero, the 
0 expression for the faulty output variance is greatly simplified, and 
Equation (6) reduces to 
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/-1 /*y/=i/=i ' ' 



Thus, under the above assumptions, hard defects do not 
modify the fault-free circuit standard deviation (3o f = 3o ;/ ) but affect only 
the mean value. 

5 Now that the fault-free and faulty circuit distributions of 

the output parameter are obtained, the detectability threshold is set to be 
the minimal distance between the fault-free and faulty circuits that 
guarantee detectability of the fault (Figure 6). This minimal distance is 
given by the relation: 

10 

D t h res h oH, = M/~M// * 30, + 3o„ (8) 



where u /; and \i f are the mean output parameter values for the fault-free 
and faulty circuit distributions, respectively, and a ff and 9 are the 
estimated fault-free and faulty output standard deviations of these 
15 distributions, respectively. 

Step 27: 

The fault value R fi could now be obtained from 
20 Equations (2) and (5): 
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R 



_ . l 3 V 3o //l 

G ou» _ our (9) 



eo ff 

R t> * ~Zo7t (10) 



combining (3) and (10), we obtain 



R A 



(11) 



out 



5 where S and G (steps 22 and 26, respectively) are obtained using the 
well known adjoint network method [Stephan W. Director and Ronald A. 
Rohrer, "Automated Network Design - The Frequency Domain Case", 
IEEE Trans, on Circuit Theory, CT-16, no. 3, August 1969, pp 330-337; 
J. Chojcan and J. Izydorcky, "The Time domain Sensitivity computation 
10 Using SPICE2. The Linear Network Case", Proceedings Intern, AMSE 
conference "Signal & System", Cetinje (Yugoslavia), Sep. 3-5, 1990, Vol 
3, pp 1 13-123; Leon O. Chua and Pen-Min Lin, "Computer-Aided Analysis 
of Electronic Circuits", Prentice-Hall, INC. "Englewood Cliffs, New Jersey, 
1975; Stephan W. Director and Ronald A. Rohrer, "The Generalized 
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Adjoint Network and Network Sensitivities", IEEE Trans, on Circuit and 
Theory, Vol. CT-16, no. 3, August 1969, pp. 318-323; and Tellegen 
B.D.H., "A General Network Theorem, with Application", Phillips Res. 
Dept. no. 7, pp 259-269]. The adjoint network method allows to compute 
the sensitivities of one output parameter with respect to all component 
5 variations (existing and non-existing components) in only two simulations, 
one for the original network and one for the corresponding adjoint 
network. The adjoint network method for sensitivity computation in AC, 
DC and transient domain has been implemented using Hspice [Meta- 
Software, "HSPICE User's Manual version H92", Meta-Software, Inc. 

10 1 992] as a basic simulator [Nairn B-Hamida, Khaled Saab, David Marche 
and Bozena Kaminska, "FaultMaxx: A Perturbation Based Fault Modeling 
and Simulation for Mixed-Signal Circuits", IEEE Asien Test Conference, 
October 1 997; Saab Khaled, Nairn B. Hamida, David Marche and Bozena 
Kaminska, "LIMSoft: Automated Tool for Sensitivity Analysis and Test 

1 5 Vector Generation", IEEE Proceedings on Circuits, Devices and Systems, 
December 1 996; and Nairn B. Hamida, Khaled, Saab, David Marche, B. 
Kaminska and Gay Quesnel, "LIMSoft Automated Tool for Design and 
Test Integration of Analog Circuits", International Test Conference, 1996]. 

20 Knowing the component tolerance AXj, the output 

parameter distribution of the fault-free and faulty circuit 307, and 3o, are 
estimated. Then equation (1 1) enables estimation of the resistance value 
R feu8 that will cause the output parameter to go out of the tolerance range. 

25 Step 28: 

The values R feu# are stored in the fault dictionary. 
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In summary, from the fault-free circuit distribution, the 
mean value and standard deviation of the output parameter due to circuit 
component variations is computed. Then, from the fault-free circuit 
distribution of the output parameter and the gradient values calculated in 
step 26, the resistance value that will drive the output parameter out of its 
5 tolerance range is computed and stored. Note that the resistance value 
which is obtained indicates the value of the resistance which, if added to 
the branch or circuit, will cause the output parameter to go out of its 
tolerance range. This resistance value will be used for test vector 
specification. 

0 

The following section describes an algorithm which uses 
the fault dictionary generation approach and fault dominance concept to 
derive the fault coverage and the test vector specification that detect the 
largest set of faults without targeting individual faults. 

5 

Fault Dominance: 

In digital circuits, fault dominance is used to reduce the 
number of faults that need to be considered. 

Definition [M. Abramovici, MA Breuer, A.D. Friedman, 
"Digital Systems Testing and Testable Design", IEEE press, 1990.]: Let 
T g be the set of all tests that detects a fault g. Let's say that a fault / 
dominates a fault g if / and g are equivalent under T r 

5 

In other words, if / dominates g then any test that 
detects g , will also detect / (on the same primary output). Therefore, for 



WO 08/55880 PCT/CA98/00538 



21 



fault detection it is unnecessary to consider the dominating fault /, since 
by deriving a test to detect g we automatically obtain a test that detect / 
as well. 

In analog circuit, the input/output relationship is more 
5 complex, but to the first order approximation, the above fault dominance 
theorem could be used as well. 

Indeed, instead of testing for the upper and lower 
resistance values for faults as described in R.J .A. Harvey, A.M.D. 
0 Richardson, E.M.J.G. Bruls, K. Baker, "Analog Fault Simulation Based on 
Layout Dependent Fault Models", ITC 95, pp-641-649 (Table 1), in the 
present invention, fault dominance will be used where only the least 
dominating upper resistance will be tested for shorts, and the least 
dominating lower resistance will be tested for opens. 

5 

Test Vector Specification: 

An example of hard fault test vector specification 
algorithm is illustrated in Figure 8. 

0 

Step 31 : A set of stimuli is selected. A default stimulus could be 

used or it can be made by the test engineers in an 
interactive mode in order to consider any special 
characteristic of a circuit. Stimuli are divided into DC, 
AC and transient stimuli [Majoj Sachdev, U A Realistic 
Defect Oriented Testability Methodology for Analog 
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Circuits", JETTA 1993.]: sine wave, pulse, ramp, any 
PWL function, etc. 

Step 32: A stimulus T is selected from the set of stimuli of step 

31- 

5 

Step 33: A fault is selected from a fault list 48 constructed by the 

algorithm of Figure 7. 

Step 34: Step 34 is responsive to the selected stimulus from step 

10 32 and the fault selected in step 33 to obtain a value of 

Rfauit fr om a fault dictionary 49. 

Step 35: Step 35 determines whether the fault is an open or a 

short. 

15 

Steps 36-40: - If the fault is an open (step 35) and R tautt < R^, (step 
36; Table 1 presents a list of circuit defects and the 
corresponding typical resistance range obtained by 
statistical analysis of circuit defects [R.J.A. Harvey, 

20 A.M.D. Richardson, E.M.J.G. Bruls, K. Baker, "Analog 

Fault Simulation Based on Layout Dependent Fault 
Models", ITC 95, pp-641 -649], the fault is removed from 
the fault list (step 39) and the fault is marked as being 
detected by simulus T (step 40). Then, it is concluded 

25 that Rfc Ult (or higher resistance value) will cause the 

output parameter to deviate out of tolerance and the 
simulus T is accepted as a valid test vector 
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- if the fault is an open (step 35) and R (auIt > (step 

36) , the fault is marked as being undetected by simulus 
T (step 38). Then, it is concluded that no deviation of 
the output parameter could be detected and the 
stimulus T is rejected as a valid test vector; 

- if the fault is a short (step 35) and R feuIl > (step 

37) , the fault is removed from the fault list (step 39) and 
the fault is marked as being detected by simulus T (step 
40). Then, it is concluded that (or higher resistance 
value) will cause the output parameter to deviate out of 
tolerance and the simulus T is accepted as a valid test 
vector; and 

- if the fault is a short (step 35) and R fautt < R^, (step 
37), the fault is marked as being undetected by simulus 
T (step 38). Then, it is concluded that no deviation of 
the output parameter could be detected and the 
stimulus T is rejected as a valid test vector. 

In this step, the algorithm is returned to step 33 until all 
faults in the fault list have not been processed for a 
given stimulus T. 

In this step, the algorithm is returned to step 32 as long 
as all the stimuli of the set selected in step 31 have not 
been processed. 
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Step 43: The fault coverage is printed; this is a list of faults that 

can be detected by the set of stimuli selected in step 
31. 

Step 44: The test vectors are printed; this constitutes the test 

5 ~ vector specification. 



Step 45: The test algorithm is terminated. 

Therefore, the fault values R faull corresponding to all the 
10 stimuli and all the faults are compared with typical values for shorts and 
opens (Table 1). More specifically, all the stimuli are evaluated for each 
fault; this iteration loop enables specification of test vectors that 
maximizes the observability of the fault. 

15 Experimental results: 

The method according to the invention were applied to 
a second order band-pass filter as shown in Figure 9, to the fifth order 
chebychev band-pass filter of Figure 12, and to several other circuits 
20 listed in Table 4. 

Example 1 (second order band-pass filter of Figure 9): 

In the filter of Figure 9, Rg=220kO, R1=10O, R2=10kQ, 
25 R3=10kO, R4=10kQ, Rd=200kO, C1=1.59nf, and C2=1.59nf. 
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The gradient method has been used for the band-pass 
filter testing at the circuit's nominal frequency response, 10 kHz. For 
simplicity, three possible opens and two possible shorts circuits were 
considered: open Rg, open R2, open R3, short 3&8 and short 6&4. The 
output voltage was used as a detection mechanism (output parameter) 
and the detection threshold was set to 5% of the nominal voltage value. 
In other words, the method of the invention was used to estimate the 
resistance value R toM beyond which the signature of the output parameter 
will be modified by more than 5% (Figure 10). 

In this example, four out of the five faults have been 
detected. R^ for open Rg, open R2, and open R3 was lower than Rt^ 
and the test vector (10KHz sine wave) was marked as a valid test vector. 
For short 3&8 R iaM was higher than R. ypica , and the test vector was also 
marked as a valid test vector for this fault. However, for short 6&4, R faul , 
was lower than R^^ and the test vector was rejected as a non-valid test 
vector for this fault. 

To validate the efficiency of this method, each fault was 
reinjected manually into the circuit and simulated. Figure 11 shows the 
results of the fault simulation over the considered catastrophic defects. 
In this figure, each fault is analysed to confirm whether or not it was 
detected at the output node at the 10KHz test frequency. Indeed, all the 
accepted resistances modified the signature of the fault-free output 
voltage by more than 5% which has been determined as a detection 
criteria. 

Example 2: Fifth order chebychev filter of Figure 1 2 



WO 98/55880 A 4* PCT/CA98/00538 



26 



The fifth order chebychev filter of Figure 12 has been 
tested for 25 possible hard faults. The set of stimuli T was selected as 
the frequency range 0-20kHz. The detection threshold was set to 5% of 
the nominal output voltage, R^, for open set to 1MQ, and ^ for 
shorts set to 1kO. 

5 

It was found that all the relevant possible shorts and 
opens are easily detected. The following Table 2 and Table 3 give some 
possible hard faults, their computed nominal values and the frequency at 
which they are detected. 

10 

Table 2 Fifth order chebychev filter open circuit fault results 



OPEN R lauM (ohm) R^ pieal (ohm) Frequency(Hz) Decision 



R1 


1.1K 


1Meg 


300 


Accepted 


R2 


1.93K 


1Meg 


300 


Accepted 


R3 


400 


1Meg 


0 


Accepted 


R4 


662.5 


1Meg 


0 


Accepted 


R5 


4.83k 


1Meg 


100 


Accepted 


R6 


2.78k 


1Meg 


200 


Accepted 


R7 


2.78k 


1Meg 


200 


Accepted 


R8 


431 


1Meg 


0 


Accepted 


R9 


555 


1Meg 


0 


Accepted 


C1 


5.12K 


1Meg 


100 


Accepted 


C2 


4.43k 


1Meg 


100 


Accepted 


C3 


2.89k 


1Meg 


100 


Accepted 


gm1 


872 


1Meg 


200 


Accepted 


gm2 


243 


1Meg 


0 


Accepted 


gm3 


1k 


1Meg 


1k 


Accepted 



Table 3 Fifth order chebychev filter short circuit fault results 



WO 98/55880 




PCT/CA98/00538 



27 



SHORT Rt,u»(ohm) R^^ohm) Frequency(Hz) Decision 



Vtn&O 


INF 


1k 


all frequencies 


Accepted 


VDD&O 


INF 


1k 


all frequencies 


Accepted 


VSS&O 


INF 


1k 


all frequencies 


Accepted 


R1&0 


20 Meg 


1k. 


0 


Accepted 


R1 &R2 


8 Meg 


1k 


0 


Accepted 


R5&R6 


85 k 


1k 


0 


Accepted 


R9&0 


50 M 


1k 


1k 


Accepted 


in- 8.0 


25 Meg 


1k 


0 


Accepted 


in+8.0 


2e-5 


1k 


0 


Rejected 



15 The obtained fault coverage is 96% with only 5 test 

vectors (0Hz, 100Hz, 200Hz, 300Hz, 1KHz). 

Example 3: 

20 A set of 7 benchmark circuits were simulated. The 

benchmark circuits ranged from a simple operational amplifier to a 
complex 8 bit current DAC (digital-to-analog converter). Level 3, Level 28 
and level 49 transistor MOSFET models were used in the sensitivity 
computation environment. 

25 

In all experiments the output voltage measurements 
and/or output current measurements were considered for testing. The 
detection threshold was set to 6o where o is the standard deviation of the 
fault-free circuit. Table 4 describes the circuit type, the test domain, the 
30 number of transistors in the circuit, the number of faults in the fault list 
obtained through schematic-based fault dictionary, and the total CPU time 
on a SPARC 10 workstation. The approximated CPU time is given for 
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serial methods using one simulation per fault. The obtained fault 
coverage is also provided for comparison purposes. 



Table 4 - Benchmark results 



10 



15 



Circuit 


Test 
domain 


Number of 
transistors 


Number 
of faults 

f 


TestMaxx Serial method Fault 
CPU time(s) CPU time(s)* coverage 


Inverter 


DC 


9 


47 


0.28 


5.68 


91.5% 


low pass filter 


Transient 


9 


51 


0.43 


7.92 


94.1% 


State variable filte 


AC 


36 


184 


3.64 


334.88 


85.3% 


Chebychev filter 


AD 


27 


149 


5.32 


348.16 


87.2% 


4 bit ADC 


DC 


153 


737 


5.46 


1341.32 


85.7% 


8 bit ADC flash 


DC 


63 


295 


7.83 


256.65 


87.4% 


8 bit current DAC 


Transient 


132 


669 


6.81 


1518.63 


98.6% 


'estimated based on one simulation per fault 











20 



The average fault coverage was 90% with small CPU 
cost. There was no way to compare the fault coverage and the simulation 
time with other publications due to the large variety of analog 
benchmarks. 



25 



Figure 1 3 presents the fault coverage as function of the 
number of test vectors, for an inverter, a low pass filter, a state variable 
filter, a chebychev filter, a 4-bit analog-to-digital converter, an 8-bit 
analog-to-digital converter flash, and an 8-bit current digital-to-analog 
converter. 
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To summarize, an automatic tool for parallel fault 
simulation and test vector generation based on a cause-effect approach 
has been presented. The concept, is new, and from our simulation 



v 



98/55880 



PCT/CA98/00538 



29 



results, we can conclude that this method is highly efficient. Indeed, fault 
coverage was as high as 98.6% with simulating time several order of 
magnitude less than the serial methods. 

A plurality of articles have been referred to in the present 
- disclosure. — Thercontents and subject matter of these articles is fully 
incorporated herein by reference. 

Although the present invention has been described 
hereinabove by way of a preferred embodiment thereof, this embodiment 
can be modified at will, within the scope of the appended claims, without 
departing from the spirit and nature of the subject invention. 
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WHAT IS CLAIMED IS: 

1. A method for constructing a fault dictionary for a 
circuit in view of subsequently testing the circuit, comprising the steps of: 

( a) describing the circuit; 

( b) extracting a list of faults from the description of the 

circuit; 

( c) in response to the description of the circuit, 
calculating a fault-free circuit distribution of an output parameter of the 
circuit; 

( d) calculating a faulty circuit distribution of the output 
parameter in response to the faults of said list; 

( e) calculating a fault value from the fault-free and faulty 
circuit distributions; 

( f) storing the calculated fault value in the fault dictionary 
in view of subsequently specifying at least one test vector for application 
to the circuit in view of testing said circuit. 

2. A method for constructing a fault dictionary as recited 
in claim 1, in which step ( b) comprises extracting from the description of 
the circuit layout-based faults to form said list 

3. A method for constructing a fault dictionary as recited 
in claim 1, in which step ( b) comprises extracting from the description of 
the circuit schematic-based faults to form said list. 
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4. A method for constructing a fault dictionary as recited 
in claim 1, wherein the circuit comprises a plurality of components and 
wherein step ( c) comprises: 

calculating sensitivities of the output parameter to 
variation of the components of the circuit; and 

calculating the fault-free circuit distribution in relation to 
the calculated sensitivities. 

5. A method for constructing a fault dictionary as recited 
in claim 1, wherein step ( d) comprises the steps of: 

selecting a stimulus; and 

calculating gradients of the output parameter with 
respect to the faults for the selected stimulus. 

6. A method for constructing a fault dictionary as recited 
in claim 5, further comprising the step of repeating steps ( d), ( e) and (f) 
for a plurality of stimuli. 

7. A method for constructing a fault dictionary as recited 
in claim 5, in which step ( e) comprises calculating the fault value in 
relation to the fault-free circuit distribution and the calculated gradients ; 

8. A method for constructing a fault dictionary as recited 
in claim 1, wherein the circuit comprises a plurality of components, and 
wherein: 
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step ( c) comprises calculating sensitivities of the output 
parameter to variation of the components of the circuit; 

step ( d) comprises calculating gradients of the output 
parameter with respect to the faults; and 

step ( e) comprises calculating the fault value in relation 
5 to the calculated sensitivities and gradients. 

9. A method for constructing a fault dictionary as recited 
in claim 1 , wherein the fault value is a minimal fault value which, if added 
to the circuit, causes the output parameter to go out of a predetermined 

10 tolerance range. 

10. A method for constructing a fault dictionary as 
recited in claim 8, wherein the fault value is a minimal resistance value 
which, if added to the circuit, causes the output parameter to go out of a 

15 predetermined tolerance range. 

11. A method for constructing a fault dictionary as 
recited in claim 10, wherein the minimal resistance value is given by the 
relation: 

20 
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wher © - SJ* is the sensitivity of the output parameter to 

variations of the circuit component x,; 

- is a standard deviation of the component 

- S*, 0 "* is the sensitivity of the output parameter to 
variations of the component x,; 

5 - Oxjj is a covariance term; and 

- Gf, 0 "' is the gradient of the output parameter with 
respect to the fault f,. 

12. A method for constructing a fault dictionary as 
10 recited in claim 1 , further comprising the steps of: 

( g) comparing the fault value calculated in step ( e) with 
a typical fault value; 

( h) for each stimulus determining whether each fault is 
detected or undetected by said stimulus in accordance with the result of 

15 said comparison; 

( i) building test vectors from the stimuli of said set most 
susceptible to detect faults in said circuit. 

13. A method for building test vectors in view of testing 
20 a given circuit, comprising the steps of: 

( a) selecting a set of stimuli; 

( b) selecting a stimulus of said set; 

( c) selecting a fault from a list of faults of said circuit; 
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( d) from a fault dictionary, finding a fault value related 
to the selected fault and stimulus; 

( e) comparing the faut value found in step ( d) to a 
typical fault value; 

( f) determining whether the selected fault is detected or 
5 undetected by applying said stimulus to the circuit including the selected 
fault in accordance with the result of the comparison of step ( e); 

( g) repeating steps ( b), ( c), ( d), ( e) and ( f) for each 
pair of stimulus and fault; and 

( h) building test vectors from the stimuli of said set most 
1 0 susceptible to detect faults in said circuit. 

14. A method for building test vectors as defined in 
claim 13, in which, when the fault is an open circuit, step ( f) comprises 
determining that the selected fault is detected by the selected stimulus 

1 5 when the result of the comparison of step ( e) is that the fault value found 
in step ( d) is lower than the typical fault value. 

15. A method for building test vectors as defined in 
claim 13, in which, when the fault is a short circuit, step ( f) comprises 

20 determining that the selected fault is detected by the selected stimulus 
when the result of the comparison of step ( e) is that the fault value found 
in step ( d) is higher than the typical fault value. 
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Intermediate list J 

J open net(T) poly defect: ((2210.0 355.0) (2390.0 535.0)) 
short netfY A") poly/poly-contact defect: ((2060.0 505.0) 
(2240.0 685.0)) 

open net(Y) poly defect: ((2210.0 505.0) (2390.0- 685.0)) 
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elumination 
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i open netfA] Weight: 10.000000 i 
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